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Abstract 

Theoretical predictions for inclusive semileptonic B decay rates are rewritten 
in terms of the T(15) meson mass instead of the b quark mass, using a mod- 
ified perturbation expansion. This method gives theoretically consistent and 
phenomenologically useful results. Perturbation theory is well behaved, and 
the largest theoretical error in the predictions coming from the uncertainty in 
the quark mass is eliminated. The results are applied to the determination of 
\Vd,\, \V ub \, and A x . 
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Inclusive decay rates of hadrons containing a heavy quark can be systematically expanded 
in powers of a s (mg) and Aqcd/^q, where mq is the mass of the heavy quark and Aqcd 
is the nonperturbative scale parameter of the strong interactions. In the uiq — > oo limit, 
inclusive decay rates are given by free quark decay and the order Aqcd/^q corrections 
vanish The leading nonperturbative corrections of order AQ CD /mg are parameterized 
by two hadronic matrix elements These results are now used to determine the CKM 

matrix elements \V c b\ and using experimental data on inclusive semileptonic B meson 
decays. 

At present, the largest theoretical uncertainties in the B — ► X c eu and B —>■ X u ev decay 
rates arise from poor knowledge of the b quark mass. The b quark pole mass is an infrared 
sensitive quantity which is not well defined beyond perturbation theory ||. This is related 
to the bad behavior of perturbative corrections to the inclusive decay rate when it is written 
in terms of the pole mass f|0]. The decay rate has been rewritten, with the hope of reducing 
the theoretical uncertainties, in terms of other quantities such as the B meson mass and the 
A parameter of HQET, or in terms of the infrared safe MS mass of the b quark. Nonetheless, 
the uncertainties remain sizable and are a significant part of the present theoretical errors 
on \V cb \ and \V ub \. 

In this paper, the theoretical predictions for semileptonic B decay rates are rewritten in 
terms of the T(15) meson mass rather than the b quark mass. This eliminates the uncertainty 
due to the ml factor in the decay rates, and at the same time improves the behavior of the 
perturbation series. Our formulae relate measurable quantities to one another and the 
resulting perturbation series is free of renormalon ambiguities. 

The inclusive decay rate B — > X u ev is 
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Here m& is the b quark pole mass, (3q = 11 — 2n//3 is the first coefficient of the QCD (3- 
function, and a s is the running coupling constant in the MS scheme at the scale ft = m b . 
The variable e = 1 denotes the order in our modified expansion. There is a subtlety in the 
power counting for the T mass, for which the difference between powers of a s and e will be 
important. Only the part of the a^' 3 corrections proportional to (3q 2 (the BLM piece M) is 
known. It is the dominant part of the two-loop correction in examples where the entire two- 
loop correction has been computed (see, e.g., Eq. (||)). The 1/m 2 terms are a few percent, 
so the (Xs/m 2 and 1/m 3 corrections are negligible. With a s (mb) = 0.22 and rif = 4, the 
perturbative series in Eq. (D) is 1 — 0.17e — 0.13blm£ 2 — 0.12 B Ltvie 3 — • • where we have 
used the subscript BLM to remind the reader that only the BLM piece of the a^' 3 terms 
has been computed. It is difficult to estimate T(B — > X u ev) reliably, since uncertainties in 
ml and in the perturbative expansion are large. Moreover, the perturbation series at large 
orders contains a contribution of order a"/?Q ~ l n\ which is not Borel summable, leading to a 
renormalon ambiguity. 

The pole mass m^ is an infrared sensitive quantity. It can be related to an infrared safe 
mass such as the MS mass via (for rif = 4) || 
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1 + ^6 + (1.56/% - 1-07)4 + • • • ( 2 ) 



This relation also has terms of the form a™/3o _1 n! at high orders. There is a cancellation 
between the a"/3o _1 n! terms in Eqs. (|l|) and (0) when the inclusive decay rate is rewritten in 
terms of the MS mass [10] . While this cancellation is present at high orders, the perturbation 



series in Eq. (JJ) with m h — > m b is 1 + 0.30e + 0.19BLMe 2 + 0.05BLMe 3 0, so there are still 
large corrections at low orders. Furthermore, using the MS mass does not remove the quark 
mass uncertainty in the decay rate. 

A simple method of avoiding problems with the quark mass is to use instead the hadron 
mass. Unfortunately, the B meson and h quark masses differ by order Aqcd ; and so this 
reintroduces a AQ CD /m 6 correction to the inclusive decay rate. A better method is to rewrite 
expressions like Eq. (fl|) in terms of the T mass to obtain well defined formulae for B decay 
rates in terms of my. The resulting expressions are free of renormalon ambiguities, and they 
express one measurable quantity in terms of another. We will also see numerically that the 
a s corrections are small when the B decay rate is written in terms of the T mass. 

There is an interesting theoretical subtlety in the behavior of the perturbation series for 
the T mass in terms of the quark pole mass. This is simplest to illustrate in the large /3q 
(i.e., bubble summation) approximation. Schematically, the perturbative expansion of the 
T mass in terms of nib is 

2^-1- {as ° F)2 h + ^U + l) + (^) 2 (f +1+1) + .. . 

+ \t + t~ 1 + . .. + !)+. ..}, (3) 

where t = ln[/i/(mf,a s CV)], Cf = 4/3, and the precise coefficients are not shown. At low 
orders this series is of the form {a 2 , alfio, ck^q, . . .}, whereas the corrections in Eqs. ([!]) 
and (0) are of order {a s , ct^A), o? s (3q, . . .}. An explicit calculation using the Borel transform 
of the static quark potential [|TTJ shows that this mismatch disappears at higher orders. The 



terms in Eq. (|]) of the form {t n + £ n ~~ 1 + . . . + 1) exponentiate to give exp(£) = /i/ (mbasCp), 
and correct the mismatch between the powers of a s and (3q. This has to happen since my is 
a physical quantity, so the renormalon ambiguities must cancel in Eq. (|3]) between 2mb and 



the potential plus kinetic energies [12 



The expression for the T mass in terms of is 
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+ (^) 2 (3. + 9f + 2 C,3,4 + IIy + ... 

The ellipsis denote terms of order a A s with at most one power of (3q or f3\ (which are known), 
as well as terms of order a 5 s . The arguments following Eq. @ show that to ensure the 
cancellation of renormalon ambiguities when we combine Eqs. ([I]) and (^), terms of order 
a™ in Eq. (^) should be viewed as if they were only of order a™ -1 . For this reason, the 
power of e in Eq. (HI) is one less than the power of a s . One should also choose the same 
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renormalization scale, /i, in Eqs. ([I]) and (|j). With this prescription, it is also expected 
that the infrared sensitivity present separately in Eqs. (H) and @ will cancel to all orders 
in perturbation theory in e. For /i of order m^, Eq. (^) shows no sign of convergence; for 
H = m b it yields m T = 2m 6 (l - O.Olle - 0.016e 2 - 0.024 BL Me 3 -...)■ Th e bad behavior of 
this series is unimportant, since the only physical question is what happens when we use 
Eq. (f|) to predict B decay rates in terms of mr- 

An important theoretical uncertainty in applying the above approach is the size of non- 
perturbative corrections to Eq. (|]). The dynamics of the T system can be described using 
NRQCD [T3]. The leading nonperturbative corrections to my arise from matrix elements in 
the T of -flight; the Hamiltonian of the light degrees of freedom. In B mesons, the leading 
nonperturbative correction to the B meson mass is due to the matrix element of ifiight, 
which is the A parameter of order Aqcd- The Aqcd dependence is different for the T. 
/flight is the integral of a local Hamiltonian density, i/ught = Jd 3 x Tiu g ht (x) . The radius of 
the T is a ~ l/(m(,a; s ), so the matrix element of i/u g ht is of order a 3 AQ CD , by dimensional 
analysis. (Note that the matrix element of Taught is of order Aq CD , not m\. Terms that 
grow with rrib can be treated using NRQCD perturbation theory.) Using l/a ~ 1 GeV, and 
Aqcd ~ 350 MeV, of order a constituent quark mass, gives a nonperturbative correction of 
15 MeV. Using instead A QCD ~ 500 MeV gives a correction of 60 MeV. We will use 100 MeV 
as a conservative estimate of the nonperturbative contribution to my. 

Substituting Eq. (§]) into Eq. ([J) and collecting terms of a given order in e gives the 
B — > X u ev decay rate in the large /3q approximation in terms of the T mass, 
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using \i = rrib and a s (rrib) = 0.22. The non-BLM parts of the e 2,3 terms have been neglected. 
The perturbation series, 1 — 0.115e — 0.035blmc 2 — 0.005BLMe 3 , is far better behaved than 
the series in Eq. (|1|), 1 — 0.17e — 0.13blmc 2 — 0.12blmc 3 , or the series expressed in terms of 
the MS mass, 1 + 0.30e + 0.19blm£ 2 + COSblmc 3 - The uncertainty in the B decay rate using 
Eq. (|5p is much smaller than that in Eq. (JIJ), both because the perturbation series is better 
behaved, and because the T mass is better known (and better defined) than the b quark 
mass. 

The non-BLM order a 2 corrections to b decay have only been calculated for b — > c decay, 



at three values of the invariant mass of the lepton pair |L5|. Extrapolating these results 
to m c — > gives the estimate that the complete a 2 correction to b — > u decay is about 
(90 ± 10)% of the order a 2 /3o result ||. With this estimate, and including the entire e 2 term 
in Eq. (Q) gives at order e 2 

T(B -> X u ev) = (^) 5 [l - 0.1156 - (0.045 ± 0.013) e 2 

- (0.20A 2 - 0.02Ai)/GeV 2 ] , (6) 

where the error on the e 2 term is due to the ±10% uncertainty in the a 2 term in b — > u 
decay. Eq. (@) yields a relation between \Vy^\ and the total semileptonic B — > X u eu decay 
rate with very small uncertainty, 
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\V ub \ = (3.06 ± 0.08 ± 0.08) x 10 



-3 



( B(B^X u ev) 1.6ps\ 1/2 

where we have used A2 = 0.12GeV 2 and Ai = (—0.25 ± 0.25) GeV 2 . The first error is 
obtained by assigning an uncertainty in Eq. equal to the value of the e 2 term and the 
second is from assuming a 100 MeV uncertainty in Eq. (f4|). The scale dependence of \V u b\ 
due to varying /1 in the range m^jl < fi < 2rrib is less than 1%. The uncertainty in Ai 
makes a negligible contribution to the total error. It is unlikely that B(B — > X u ev) will 
be measured without significant experimental cuts, for example, on the hadronic invariant 
IB|. Our method should reduce the uncertainties in such analyses as well. 



mass 



The B — > X c ev decay depends on both m& and m c . It is convenient to express the decay 
rate in terms of mx and Ai instead of and m c , using Eq. ([|) and 

_ — ( Ai Ai \ 

m b - m c = m B - m D + — —7- + . . . , (8) 

\2m B 2m D J 

where = (3m#. + m#)/4 = 5.313 GeV and run = (3m£>* + m/))/4 = 1.973 GeV. The ot s 

can 



correction to free quark decay is known analytically ||17|| , and the full order a 2 s result |15 
be estimated numerically (at the scale /i = m h ) by multiplying the order o; 2 /?o correction 
by 0.9 ± 0.05. We then find 
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where the phase space has also been expanded in e. For comparison, the perturbation 
series in this relation when written in terms of the pole mass is 1 — 0.12e — 0.06e 2 — . . . 0. 
Equation (0) implies 



\V cb \ = (41.6 ± 0.8 ± 0.7 ± 0.5) x 10 



-3 



(B(B X c eu) 1.6ps\ 1/2 . , 
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where ?]qed ~ 1.007 is the electromagnetic radiative correction. The uncertainties come 
from assuming an error in Eq. (P) equal to the e 2 term, the 0.25 GeV 2 error in Ai, and 
a 100 MeV error in Eq. (|j), respectively. The second uncertainty is reduced to ±0.3 by 
extracting Ai from the electron spectrum in B — > X c eu; see Eq. (|Tl~D. The agreement of 
\V c b\ with other determinations (such as exclusive decays) is a check that nonperturbative 
corrections to Eq. @ are indeed small. 

In Ref. [18|] A and Ai were extracted from the lepton spectrum in B — ► X c ev decay. 



With our approach, there is no dependence on A, so we can determine Ai directly with small 
uncertainty. Considering the observable R\ = j 1 5GeV Ei(dT/dE e )dEi/ J 15GeV (dT/dE e )dE£, 
a fit to the same data yields 

Ai = (-0.27 ± 0.10 ± 0.04) GeV 2 . (11) 
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The central value includes corrections of order a 2 /?o fL9 |. The first error is dominated by 
1/m 3 , corrections [pOR . We varied the dimension-six matrix elements between ±(0.5 GeV) 3 , 
and combined their coefficients in quadrature in the error estimate. The second error is from 
assuming a 100 MeV uncertainty in Eq. (f|). The central value of Ai at tree level or at order 
a s is within 0.03 GeV 2 of the one in Eq. flTT|). 

The above results can also be applied to D —> Xev decay, using a s {m c ) = 0.35 and 
rif — 3. Nonperturbative effects are clearly much larger in the J/ijj than in the T, so one 
might expect the entire analysis to break down completely. It is remarkable that this does 
not occur. Using raj^ = 2m c (l — 0.027e — 0.059e 2 — 0.130e 3 — ...), neglecting m s , and 
following the same procedure as for b — > u decay, we find 
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The e 3 contribution to Eq. ( |T2"D is larger than the order e 2 term. The perturbation series 
expressed in terms of the pole mass has a much worse behavior, roughly 1 — 0.27e — 0.32e 2 . 
Using A 2 (m c ) = 0.14 GeV 2 and Ai from Eq. (|TT|), we obtain 

\V CS \ 2 + \V cd \ 2 = (1.00 ± 0.06 ± 0.04) 

(BjD^Xev) 1.06 p 3 \ 

V 0.17 t d ± )' v ' 

where the uncertainties come from assuming an error in Eq. ( |T2"D equal to the e 2 term and the 
error in Ai, respectively. We have not included an estimate of nonperturbative corrections 
to the J/ip mass, or of scale dependence. The LEP measurements of the hadronic W width 
yield \ V CS \ = 0.98 ±0.07 ±0.04 [^TJ. The uncertainty in Eq. (|T3| ) is comparable to this, since 
the experimental error of E{D ± — > Xev) is about 10%. Eq. ( |T3"D has theoretical uncertainties 
which we cannot estimate. The validity of quark-hadron duality may be questionable since 
the final states are almost saturated by K and K* . In addition, an estimate similar to 
that for the T suggests that the nonperturbative contribution to the J/ip mass is of order 
500 MeV (using 1/a ~ 0.5 GeV and Aqcd ~ 500 MeV). This gives an uncertainty of order 
100% in |V C;S | 2 ± |V^| 2 . The agreement of Eq. (p~3|) with the experimental results may be a 
coincidence, or may signal that nonperturbative corrections in the mass relation are much 
smaller than naive expectations. 

We have chosen to write our B decay results in terms of the Y(IS') mass. One could 
equally well write them in terms of the mass of excited states, such as the T(2S). The 
perturbation series is expected to be worse behaved than for the T(IS'). The main difference 
is in the estimate of nonperturbative corrections to the T(2S) mass. The radius of the 2S 
state is about four times that of the IS, so the nonperturbative corrections, which grow as 
a 3 , are approximately 64 times larger. This implies a similar increase in the error on the 
CKM angles. Ignoring nonperturbative corrections for the moment, the analog of Eq. (|j) 
for the T(2S') evaluated at the scale /i = nib is mr(2S) = 2m&(l — 0.0027e — 0.0059e 2 — 
0.0117blm£ 3 — ...)• Numerically, the first few corrections are smaller than for the T(IS'), 
but the convergence of the series is worse. The B — > X u ev decay rate in the large (3q 
approximation in terms of the T(2S*) mass is then 
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1 - 0.155e - 0.098 BL Me 2 - 0.065 BL Me 3 -...]. (14) 

Compared to Eq. ([5]), the convergence is worse, as expected. Nevertheless, even this formula 
gives a reasonable extraction of \V u b\. The ratio of \V u b\ 2 extracted using the 2S and IS* 
masses is [Eq. (|TJ)]/[Eq. (H)] = {1.34, 1.27, 1.17, 1.08}, where the nth number is obtained 
by truncating both equations at order e n_1 , and neglecting the Ai,2 corrections. The large 
difference at "tree level", (rrir(2S)/ m r) 5 = 1.34, is reduced by the series of perturbative cor- 
rections. Expressing the B —>■ X c ev decay rate in terms of the T(25) mass, the perturbative 
corrections in Eq. (EJ) become 



r /n x -, _ G 2 F \V cb \ 2 f m r(2S) \' 



x 0.447 



0.107e- 0.046e" 



(15) 



Again, the convergence of the series becomes worse. However, the ratio of \V c b\ 2 ex- 
tracted using the 2S and IS" masses is consistent with our estimates of the uncertainties, 
[Eq. (p3|)1/[Eq. (|)] = {1.12, 1.10, 1.08}, where the nth number is obtained by truncating 
both expressions at order e n_1 . The difference between the T(2S) and T(1S) results provide 
an estimate of nonperturbative contributions to the T mass. They suggest that nonpertur- 
bative effects are smaller than the conservative estimate we have used; they are certainly 
much smaller than the naive estimate above of a 64 x 100 MeV = 6.4 GeV nonperturbative 
contribution to the T(2S) mass. 

We have shown that inclusive semileptonic B decay rates can be predicted in terms of 
the Y(IS') mass instead of the b quark mass. It is crucial to our analysis to use the modified 
expansion in e rather than the conventional expansion in powers of a s . Our formulae relate 
only physical quantities to one another. They result in smaller theoretical uncertainties 
than existing numerical predictions, and the behavior of the perturbation series is improved. 
Moreover, the uncertainties can be estimated without resorting to cumbersome arguments, 
and they can be checked using the experimental data. 

Our main results are Eqs. (JTTJ) and (|?p, which relate the total semileptonic B — > X c u ev 
decay rates to |V^,| and \V u b\. The uncertainties are below 5% at present, and it may be 
possible to reduce them further. Our determination of Ai is given in Eq. fly]). We hope that 



applications of the method introduced in this paper will prove useful — besides reducing the 
uncertainties of \V c b\ and \V u b\ - - in analyzing a large class of data emerging from present 
and future B decay experiments. Details of our method, as well as other applications, such 
as to nonleptonic and exclusive semileptonic B decays, will be discussed elsewhere [^2] . 

We thank Mark Wise for useful discussions. This work was supported in part by the 
DOE grant DOE-FG03-97ER40546 and by the NSF grant PHY-9457911. 
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